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Abstract 

We consider an online routing problem in continuous time, where 
calls have Poisson arrivals and exponential durations. The first-fit 
dynamic alternative routing algorithm sequentially selects up to d 
random two-link routes between the two endpoints of a call, via an 
intermediate node, and assigns the call to the first route with spare 
capacity on each link, if there is such a route. The balanced dy- 
namic alternative routing algorithm simultaneously selects d random 
two-link routes; and the call is accepted on a route minimising the 
maximum of the loads on its two links, provided neither of these two 
links is saturated. 

We determine the capacities needed for these algorithms to route 
calls successfully, and find that the balanced algorithm requires a much 
smaller capacity. 
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1 Introduction 



We consider here an online routing problem in continuous time, where calls 
have Poisson arrivals and exponential durations. First, let us recall the fol- 
lowing online routing problem in discrete time from [8], where calls do not 
end. There is a set V = {1, . . . ,n} of n nodes, each pair of which may wish 
to communicate. A call is an unordered pair {u, v} of distinct nodes, that is 
an edge of the complete graph Kn on V. For each of the N = (2) unordered 
pairs {u,v} of distinct nodes, there is a direct link, also denoted by {u,v}, 
with capacity Di = Di{n). The direct link is used to route a call as long 
as it has available capacity. There are also two indirect links, denoted by uv 
and vu, each with capacity D2 = D2{n). The indirect link uv may be used 
when for some w a call {u, w} finds its direct link saturated, and we seek an 
alternative route via node v. Similarly vu may be used for alternative routes 
for calls {v,w} via u. 

We are given a sequence of M calls one at a time. For each call in turn, 
we must choose a route (either a direct link or an alternative two-link route 
via an intermediate node) if this is possible, before seeing later calls. These 
routes cannot be changed later, and calls do not end. The aim is to minimise 
the number of calls that fail to be routed successfully. 

The calls are independent random variables Zi, Z2, . . . , Zm, where each 
Zj is uniformly distributed over the edges e e E{Kn). Let d he a (fixed) 
positive integer. A general dynamic alternative routing algorithm GDAR 
operates as follows. For each call e = {u, v} in turn, the call is routed on the 
direct link if possible; and otherwise nodes wi, . . . ,Wd are selected uniformly 
at random with replacement from V \ {u,v} and the call is routed via one 
of these nodes if possible, along the two corresponding indirect links. The 
first-fit dynamic alternative routing algorithm FDAR is the version when we 
always choose the first possible alternative route, if there is one. The balanced 
dynamic alternative routing algorithm BDAR is the version when we choose 
an alternative route which minimises the larger of the current loads on its 
two indirect links, if possible. 

Results for this model were first obtained in [HIH], and later strengthened 
and extended in [8]. Consider the case where M ~ cN for a constant c > 0. 
It is known that with the algorithm FDAR we need both link capacities 



that is 'with probability — > 1 as n — 00', all M calls are routed successfully. 



Di, D2 of order 




ensure that asymptotically almost surely (aas). 
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On the other hand, the balanced method BDAR succeeds with much smaller 
capacities. Specifically, there is a tight threshold value close to In Inn/ In d 
for D2 to guarantee that aas no call fails (and the precise value of Di is 
unimportant). 

Here we consider a related continuous-time network model, with the de- 
sirable additional feature that calls end. Calls arrive in a Poisson process 
with rate AA^ = A (2), where A is a positive constant. The calls are iid 
random variables Zi, Z2, . . ., where Zj is the j-th call to arrive and is uni- 
form over the edges of Kn for each j; also let Tj be the arrival time of call 
Zj. For each edge {u,v} there are two links, uv and vu, both with capacity 
D = D{n) < 00. Since in [H] the use of direct links was found to have a minor 
effect on the total capacity requirements for efficient communication, here we 
do not use direct links but instead demand that each call be routed along 
a path consisting of a pair of indirect links. If a call is for {u,v}, then we 
pick d possible intermediate nodes uniformly at random with replacement, 
as in the GDAR algorithm. The FDAR algorithm always chooses the first 
possible alternative route, if there is one. The BDAR algorithm chooses an 
alternative route which minimises the larger of the current loads on its two 
links, if possible. Call durations are unit mean exponential random vari- 
ables, independent of one another and of the arrivals and choices processes. 
Whenever a call terminates, both busy links are freed. 

All random processes considered here are assumed to be right-continuous 
(as is standard). For each edge e = {u,v} G E and node w E V \ e, let 
Xt{e, w) denote the number of calls in progress at time t which are routed 
along the path consisting of links uw and vw, that is calls between the end 
nodes u and v of e routed via w. We call Xt = {Xt{e, w) : e E E,w E V \ e) 
the load vector at time t, and let = (Z+)^*^"~^) denote the set of all possible 
load vectors. The load vectors Xt for t > form a continuous-time discrete- 
space Markov chain. By standard results, there exists a unique stationary 
distribution 11; and, whatever the distribution of the starting state Xq, the 
distribution of the load vector Xt at time t converges to H as t — 00. 

Note that if there were no capacities then the calls in progress would form 
an immigration-death process with immigration rate A (2) and death rate 1. 
Thus in equilibrium the total number of calls in progress at time t is 

stochastically at most Po(A(2)), and in particular 




(1) 
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Our main interest is in the blocking probability, that is the probabihty 
that a new call fails to find an available route. As in the discrete version 
analysed in [8], or in the models analysed in [6] and |7], see also [2l [3l [Tl] [12] . 
we observe the 'power of two choices' phenomenon; that is, with the BDAR 
algorithm for d > 2 the capacity required to ensure that most calls are routed 
successfully is much smaller than with the FDAR algorithm. Let us now state 
our two theorems. 

Theorem 11.11 shows that, when the FDAR algorithm is used, capacity 
D(n) of order , is needed in order to ensure that no call is lost in an 

^ ' In m n 

interval of length polynomial in n. The set-up is that we have capacity 
D = D{n) ~ a ^^^^^ , and consider a time interval of length . In the case 
= 2, the result is roughly that, if a > + 1 then all calls are expected 
to be successful, and ii a < K + 1 then many calls are expected to fail. For 
general d we need to consider how K compares to c? — 2, as there is a change 
of behaviour aX K = d — 2. 

We allow any initial configuration Xq (for which E||Xo||i is finite). If the 
system is in equilibrium at time then our results apply to any time interval 
of length n^: if the system is not in equilibrium at time then we need to 
assume that the starting time to of the interval is sufficiently large. 

Theorem 1.1 Let X > be fixed and let d be a fixed positive integer. Then 
there exists a constant k = n{X, d) such that the following holds. Let K > 
and a > be constants. Let D = D{n) ~ a ^^^^ ■ Let ti = ti(n) > ln(E || 
Xolli +1) + ulnn, and consider the interval [ti,ti + n^]. 

(a) Suppose that a > {K + 2)/d if K < d — 2, and a > K + 3 — d if 
K > d — 2. Whatever version of GDAR we use, the mean number of 
calls lost during the interval is o(l). 

(b) Suppose that a < {K + 2)/d and 0<e<K + 2-daifK<d-2; 

and that a < K + 3 - d and 0<e<K + 3- d- aifK>d-2. If 
we use the FDAR algorithm, then the mean number of calls lost during 
the interval is Q{n'^). 

If Xq is in equilibrium then E ||Xo||i is finite by ([1]), so the conclusions above 
apply to an interval [ti, ti+n^] for sufficiently large ti; and hence they apply 
to all intervals [t, t + n^] for t > 0. Observe that for K < d ~ 2 the critical 
value of a is {K + 2)/d, and for K > d — 2 the critical value of a is K + 3 — d; 
and neither value depends on the arrival rate A. As foreshadowed above, the 
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next result shows that the BDAR algorithm requires significantly smaller 
capacities. 

Theorem 1.2 Let X > be fixed and let d > 2 be a fixed integer. Let K > 
be a constant. Then there exist constants k, = /t(A, d) and c = c(A, d, K) > 
such that the following holds. Let ti = ti{n) > ln(E||Xo||i +1) + n\nn, and 
consider the interval + n^]. 

(a) If D{n) > lnlnn/lnc/ + c and we use the BDAR algorithm, then the 
expected number of failing calls during the interval is o(l). 

(b) If D{n) < Inlnn/lnrf — c and we use any GDAR algorithm, then aas at 

least n^+'^-d-o{i) ^^^^ ^^^^ during the interval. 

We shall give further details, for example concerning the numbers of full 
links, when we prove Theorem 11.21 

We mention that a process similar to the one defined above, but also with 
direct links, was considered in [9], and then in [1]. The first of these works 
obtained, heuristically, some preliminary results. In [1], Anagnostopoulos 
et al. find an upper bound of In Inn/ In ci + o{ln\nn/ had) for the capacity 
required by the BDAR algorithm to ensure that, in equilibrium, an arriving 
call is accepted with probability tending to 1 as n — oo. Further, they 
identify a lower bound of (^(y^hin/liniin) for the capacity needed by the 
FDAR algorithm to achieve the same effect. Here we prove sharper versions of 
these bounds, and supplement them with a lower bound on the performance 
of the BDAR algorithm and an upper bound on the performance of the FDAR 
algorithm. Further we extend them to prove upper and lower bounds on the 
performance of these algorithms over long time intervals. In comparison 
with our bounds for the FDAR algorithm are of the order Inn/ In Inn, 
not y^hinTlnlnn; this is due to the fact that we do not allow direct routing 
between pairs of nodes. 

Let us close this section by giving some further definitions and notation 
which we shall need shortly, and then a brief plan of the paper. 

Given an edge e = {u,v} G E, let Xf(e) = ^^^f,Xt{e,w) denote the 
number of calls between u and v in progress at time t. Also, given distinct 
nodes v and u, let Xt{vu) = J2w=iu v -^tiiv, w}, u), which is the load of link vu 
at time t. Given a node v, let Xt{v) = J2u=iv ■^ti'^'^) ' "which is the number 
of calls with one end v at time t. Thus ||X(||i= is the total 

number of calls at time t. We say that a link is saturated (or full) if it has 
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load equal to its capacity D. Given a node f , we let Sf{aX v) denote the set 
of saturated links vw for calls at v at time t; and let Sfi^aX v) = |iS°(at v)\, 
which is the number of saturated links vw. Similarly, given a node w, we 
let S]^ (via w) denote the set of saturated links vw for calls at some node v 
at time t; and let S'P(via w) = \S^{via w)\. Also, for each time t we let (pt 
denote the cx-field generated by {Xg : s < t), and let 0t_ denote the a-field 
generated by {Xg : s < t). 

The rest of the paper is organised as follows. Section [2] contains some 
preliminary lemmas that will be needed in our proofs. In Section [3] we es- 
tablish concentration of measure inequalities for the total number of calls as 
well as the number of saturated links. In Section H] we prove Theorem II. 
and in Section O we prove Theorem 11.21 

2 Preliminary results 

In this section we give some elementary results which will be used in our 
proofs. 

If X has the Poisson distribution with mean /z let us write X ~ Po{fi). 
For such a random variable X we shall use the bound 

Pr(X > A;) < fi''/k\ < {eii/kf 

for each positive integer k. Given a positive integer D, let 

k>D 

be the probability that a Po{fi) random variable takes value at least D. The 
following are a pair of standard concentration inequalities for a binomial or 
Poisson random variable X with mean /i: 

Pr(X - M > e/i) < exp(-^eV) (2) 

and 

Pr(X - /i < -efi) < exp(-ieV) (3) 

for < e < 1 (see for example Theorem 2.3 (c) and inequality (2.8) in [TO]). 

We shall use the following version of Talagrand's inequality, see for ex- 
ample Theorem 4.3 in [lOj. (In the notation in [10], the function h below is 
a (c^r)-configuration function.) 
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Lemma 2.1 Let X = (Xi, X2, . . .) be a finite family of independent random 
variables, where the random variable Xj takes values in a set Qj. Let Q = 

Let c and r be positive constants, and suppose that the non-negative real- 
valued function h on Q satisfies the following two conditions for each x G ^2. 

• Changing the value of a co-ordinate Xj can change the value of h{x.) by 
at most c. 

• // h{-x) = s, then there is a set of at most rs co-ordinates such that 
h{x.') > s for any x' G which agrees with x on these co-ordinates. 

Let m be a median of the random variable Z = /i(X). For each x > 

P{Z >m + x) <2exp(-——^^—-Y (4) 



and 



x^ 



P{Z <m-x) <2expi -— . (5) 

Ac rm 



The next lemma concerns hitting times of a generahsed random walk with 
'drift'. It is the 'reverse' of Lemma 7.2 in [B]|, and can be deduced from that 
result by replacing the Yi with —Yi\ we omit the details. 

Lemma 2.2 Let 0o ^ 0i ^ • • • be a filtration, and let Yi,Y2, . . . be random 
variables taking values in { — 1,0,1} such that each Yi is (pi-measurable. Let 
Eq, El, ... , be events where Ei G (pi for each i, and let E = CiiEi. For each 
t G N, let Rt = Ro + X]!=i ^i- < p < 1/3, /ei ro and ri be integers such 
that ri < tq, and let m be an integer such that pm > 2(ro — ti). Assume that 
for each i = 1, . . . ,m, 

Pr{Yi = l|(/),_i) <p on Ei_i n > ri), 

and 

Pr{Yi = -l|0i_i) > 2p on Ei.i n > ri). 

Then 

Pr {E n {Rt > ri Vt G {1, . . . , m}}\Ro = Tq) < exp . 
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We can use the last lemma to upper bound hitting times for a type of discrete- 
time 'immigration-death' process. 

Lemma 2.3 Let (po Q (f)i ■ ■ ■ , be a filtration, and letYi,Y2,... be random 

variables taking values in {—1,0, 1} such that each Yi is (pi-measurable. Let 
Eqj . . . be events where Ei G (f)i for i = 0,1, . . .. Let a,b > be constants. 
Let r and r be integers with 2a/b < r < r — 1. 

Let Ro = ro; and let Rt — Ro + Yi. Assume that for each i — 1,2, . . . 
and each y >r -\-l 

Pr(Fi = l|0i_i) <a on Ei_^ n > r); 

and 

Pr(Fi = -l|0i_i) > by on Ei_i n = y) 

for each y — r + 1, . . . ,r, and 

Pr(Yi = -l|(/>,-i) > br on Ei^i n > f). 

Let m' = \l0g2 ^-~\, and let E be the event fl^^E'i. Then 

Pr{E n{Rt>ryte{l,..., m'}}) < 2 exp (-^) . (6) 

Proof. Let k = [log2 ^] - 1, so that 2^V <f< 2'=+V. Let Tq, Ti, . . . , be 
the hitting times to cross the k intervals from f down to 2''r, from 2^r down 
to 2''~^r, and so on, ending with the interval from 2r down to r. Thus 

To = mm{t >0:Rt = 2^}, 

and for j' = 1, . . . , k, 

Tj = mm{t > T,_i : Rt = 2^-^r}. 

Consider j G {0, . . . , A;}. We may use the last lemma with p as = b 2}^~^~^r, 
ro = 2}'~i^^r (except that for j = we let ro = f), ri = 2'^~''r and m as 
ruj = [|] . Note that Pjmj > 2'^~^+^r, which is at least twice the length of 
the interval. (It may look at first sight that we are 'giving away' rather a lot 
on the 'upward' probability but this makes only a constant factor difference.) 
Hence, with T_i = 0, 

Pr{E n {T, - T,_i > m,}) < exp (-^) < exp (-^) • 
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But now 

Pr{En{Rt> rWt E {l,...,m'}}) < 

< 
< 

Hence 

Pr(En >r\/te {l,...,m}}) < 2e~^, 

by the above if e^TS < i and trivially otherwise. I 

The next lemma is Lemma 7.3 in [6], and shows that if we try to cross 
an interval against the drift we rarely succeed. 

Lemma 2.4 Let a be a positive integer. Let p and q be reals with q > 
p > and p + q < 1 . Let 0o ^ 0i ^ 02 ^ • • • &e a filtration, and let 
Yi,Y2, . . . be random variables taking values in { — 1,0, 1} such that each 
is (pi-measurable. Let Eq,Ei, . . . be events where each Ei G (pi, and let E = 
DiEi. Let Rq = and let = Yl'i=i /^'^ ^ ~ 1, 2, . . .. Assume that for 
each i = 1, . . . ,m, 

Pr(r, = <p on Ei_^ n (0 < R,_i < a - 1), 

and 

Pr(r, = >q on Ei_i n (0 < < a - 1). 

Let 

T = mi{k > 1 : i?fc G {-l,a}}. 

Then 

Pr{En{RT = a}) < {p/qf. 

We shall require another lemma, similar to Lemma 2.1 in [6J. Consider 
the n-node case, with set VL = (Z"'")^'^""^) of all load vectors. Let us say that 
a real- valued function / on f2 has bounded increase at a node v if whenever s 
and t are times with s < t, then f{xt) is at most f{xs) plus the total number 
of arrivals in the interval {s,t] for v; f has bounded increase via a node v 
if whenever s and t are times with s < t, then f{xt) is at most f{xs) plus 



Pr{E n {Tj - Tj_i >mj)) 

j=0 

e"T4 /(I — e"i4). 
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twice the total number of arrivals in the interval (s, t] routed via v as the 
intermediate node; and / has strongly bounded increase at a node v if f{xt) 
is at most f{xs) plus the maximum number of arrivals for v in the interval 
{s, t] which use any given link incident on v. Thus for example, given v E V, 
f{x) — x{v) has bounded increase at v, f{x) — \{w & V \ {v} : x{wv) > 
D}\ has bounded increase via v, and f{x) — max„gy\{^} has strongly 

bounded increase at v. 



Lemma 2.5 Let v be a node in V . Let a,T > and let a, b be non-negative 
integers. Let / : Q — > R and g,h :'R be functions. Let E E Suppose 

that, for alia eM and all times ti <t < ti+T, Pr(^En{f{Xt) < a}^ < g{a) 
and PriE n {/(X^) > a} ) < h{a). Assume further that 



(a) f has bounded increase at v and 6 — Pr(Po(A(n — 1)0") > & + 1), or 
(h) f has strongly bounded increase atv andO — (n— l)Pr(Po(A(icr) > 6+1). 
Then 



Pr 

and 
Pr 



E^{f{Xt)<aforsomete[t^M + r]}] < {- + !) {g{a + b) + 6) , (7) 

J (7 

Er\{f{Xt)>a + b for some t e[tiM + r]}\ < (- + 1) (/i(a) + ^) . (8) 

J (7 



(c) Suppose that f has bounded increase via v and 6 — Pv{Po{\d{n — 
l)a/2) > 6 + 1). Then 



Pr 

and 
Pr 



E^{f{Xt) < a for some t e [ti, ti+r]} < (^+1) {9{a + 2b) 



(9) 



Er\{f{Xt) > a + 2b for some t e [ti,ti + T]} < {- + 1) {h{a) + 9) . 

(10) 



Proof. Consider first the case (a), when / has bounded increase at v. Note 
that the j = [^\ + 1 disjoint intervals [ti + (r — l)a, ti + ra) for r = 1, . . . , j 
cover [ti,ti + t]. Let Ar denote the event that there are at least b + 1 
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arrivals for node v in the interval [ti + (r — 1)(T, ti + rcr), so that Pr(yl,.) = 
Pr[Po(A(n - > 6 + 1] = 0. But 

E n {/(Xt) < a for some t G [ti, ti + r]} 
C E n I (u^,=i{/(Xi,+(,_i).) < a + 6}) U (U^,=i^)}, 

and (171) follows. Similarly 

E n {f{Xt) > a + 6 for some t G [ti, ti + r)} 

c E n { (u^-;,{/(x,,+er_i).) > a}) u [vjU^r) }, 

and ([8]) follows. 

To handle the case (b) when / has strongly bounded increase at f , note 
that the arrival process onto any given link vu is stochastically dominated 
by a Poisson process with rate 

in - 2 / < \d. 

^ ' {n-2Y 

Thus if Br denotes the event that there are at least 6+1 arrivals in the 
interval \ti + (r — l)cr, ti + ra) that are routed on some link vu, u ^ v, then 

Pr(5^) <{n- l)Pr[Po(Ad(T) > 6 + 1]; 

and we can complete the proof as above. 

Finally, in the case (c) the arrival process onto links with v as the inter- 
mediate node is stochastically dominated by a superposition of ("2^) inde- 
pendent Poisson processes, each with rate 



{n-2Y ~ n-2 

If Cr denotes the event that there are at least 6+1 arrivals in the interval 
[ti + (r — l)cr, ti + rcr) that are routed via v, then Pr(C,.) < Pr[Po{Xd{n — 
l)o"/2) > 6 + 1]. The rest of the proof is as above. ■ 

Consider a continuous-time Markov process (X^) with countable state 
space S and with g-matrix q = {q{x,y) : x,y E S). Under certain con- 
ditions we can compare features of its behaviour with that of independent 
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immigration-death processes. We shall need the following lemma to handle 
the lower bound part of Theorem II. 1[ 

Let be a positive integer and let the index j run over {1, . . . , A^}. For 
each j let ej denote the jth unit iV-vector and let fj be a function from 5" 
to the non-negative integers; and write f{x) for {fi{x), . . . , Jn^x)). Assume 
that the following two conditions hold: 

(i) for all distinct x and y in S such that q{x, y) > Owe have f (y) = f (x) ie^ 

for some j ; and 

(ii) for each x G 5* and each j 

X] (l{x,y) = fj{x). 

y(^S:f,{y)=f,{x)~l 

Now define Xj{x) for each x G S* and each j by setting 

yeS:fjiy)=f,{x)+l 

Lemma 2.6 For each j let Xj > be a constant. Let < ti < ^2- Let 
Yt = {{Yj)t : j = 1,...,N) be a vector of independent immigration- death 
processes where {Yj)t has immigration rate Xj and death rate 1, and has 
population at time ti . Let F ^ S be such that for each x F and each j 
we have Xj{x) > Xj, and let A be the event that Xt E F for each t G [^1,^2]- 
Then for each downset B in {0, 1, . . , 

Pr({f (X^J G i?} n A) < Pr(Y,, G 5). 

Now let Hj be a given positive integer for each j. Let = {{Yj)t : j = 
1,...,N) be like Yt except that (^)t has upper population limit nj. Let 
F O S be such that, for each x E F and each j = 1, . . . ,N, if fj{x) < nj 
then Xj{x) > Xj. Let A be the event that Xt E F for each t G [ti,t2]- Then 
for each downset B in {0, 1, . . , 

Pr({f(XiJ G 5} n i) < Pr(Y,, G S). 
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Proof. Let Xq G F, and condition on Xf^ = xq. Then we may assume that 
^j{x) > \j for each x G S", since the values \j{x) for x ^ F are irrelevant; 
and then we may ignore the event A. But now an easy coupling argument 
shows that 



and since this is true for each xq E F the result follows. The second part, 



3 Failure probability, total load and saturated 



In this section we give upper and lower bounds on the failure probability of 
a call, upper and lower bounds on the total number of calls for a node v, 
and upper bounds on the number of saturated links incident on v over long 
periods of time. All the results are valid for any GDAR algorithm. 

First we consider the failure probability of a call. For each time t and 
node V, let dt{v) = Sf{at v), and let = maxy dt{v). Recall that T is 
the random departure time of the last one of the H-'^oHi initial calls, and for 
j = 1, 2, . . . the call Zj arrives at time Tj. 

Lemma 3.1 On Tj — t 



Pr(f(XiJ G B\Xt, = xo) < Pr(Y,, G B) 



with population limits, may be proved similarly. 



links 




(11) 



and 




(12) 



also, 




(13) 



and so 




(14) 
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Proof. Recall that = (2) . On Tj = t we have 



Pr(Z, fails < -)^( ^-^ J 

and both upper bounds follow. (For the second inequahty we used the fact 
that (x + yY < 2'^-\x'^ + y^) for x,y > 0.) 
On the other hand, on Tj = t 

Vv{Zt fails > l(n-2)-"5^5^K_(i;) -Ix._m=d)'. 
But for each v & V , 

5^K4^)-W(™)=i^)' = {n-l-dUv))dt^{vY+dt^{v){dt_{v)-lY > ]^{n-l)dUvf 

for n > 3. Hence on Tj = t 

Pr(Z, fails |0,_) > l(n-2)-'^5^^(n-lK_(t;)'^ 

vev 

vev 

and so both lower bounds (fT3|) and (fT4|) follow. ■ 

To obtain our estimates for the total number of calls for a node v, and 
upper bounds on the number of saturated links incident on f , we compare 
the process {Xf) to a 'superprocess' (Xt) which satisfies Xq = Xq and evolves 
as follows. The unordered pairs of distinct nodes u and v receive independent 
rate A Poisson arrival streams of calls; each link uv has infinite capacity; and 
each call throughout its duration occupies d two-link routes chosen uniformly 
at random with replacement. (If a route is chosen more than once by a given 
call, the call will still be counted only once on the corresponding two links.) 
All call durations are unit mean exponentials independent of one another 



14 



and of the arrivals and choices processes. For each pair of distinct nodes 
u and V, Xt{{u,v}) = J^w^uv-^tii'^^'^}^''^) denotes the number of calls in 
progress between u and v at time t. Note that the process {Xt{{u,v}) : 
u,v E V,u v)t>o is itself Markov, since the capacities are infinite. It has 
a unique equilibrium distribution, and in equilibrium the Xt{{u,v}) are all 
independent -Po(A) random variables. Thus, in equilibrium the total number 
||Xt||i of ongoing calls at time t is Po{X(^)); and, for each v, the total number 
Xt{v) of ongoing calls with one end v is Po{\{n — 1)). 

We shall use T„ to denote the time that the last of the Xq{v) initial calls 
with one end v departs. Also, we let T = max^gyT^. For various events A 
we shall give an upper bound on Pr(y4 fl {T < t}). We may later obtain an 
upper bound on Pr(y4) using 

Pr(A) < Pr(An {T < t}) +Pr(T > t), (15) 

and noting that 

Pr(r > t) < E||Xo||i e-*. (16) 

To see why f fT6|) holds, temporarily let St be the number of initial calls 
surviving to time t, and observe that 

Pr(T >t) = Pr(5t > 0) < E5t = E||Xo||i e"* 

The next lemma shows that for any node v E V , Xt{v) is unlikely to 
deviate far above A(n — 1) once the initial calls have gone. 

Lemma 3.2 Let < 6 < 1, let n be a positive integer, and let At be the 
event that Xt{v) > (1 + (5)A(n — 1) for some vertex v. Then for all times 
t > 

Pr{At n{T< t}) < ne-^-^'^^"-^). (17) 
(The value of D is not relevant here.) 

Proof. Consider links vu incident on a given node v. The total number of 
calls on those links is stochastically dominated by the number of calls for 
V corresponding to the process (Xt). Let (Yt) = {Yt{{u,v}) : u,v E V) for 
i > be a Markov process with the same g-matrix as (Xt) but in equilibrium. 
We couple (Xt), (Xt) and (Yt) as follows. We assume that Xq = Xq. All 
subsequent arrival and potential departure times of new calls are the same 
for the three processes, except that the departures of calls that were not 
accepted in (Xt) are ignored in that process. Additionally, every one of the 
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||Xo||i initial calls in (Xj) is coupled with a corresponding initial call in (Xf), 
and the paired calls have the same departure times. 

Under the coupling, on the event < t, for all times t and all m 7^ f , 

Xt{vu) < Xt{vu) < Yt{vu), 

and so also 

Xtiv) < Xtiv) < Ytiv). 

But Yt{v) is a Poisson random variable with mean A(n — 1), and so by the 
concentration inequality ([2]), we have 



Pr {{Xt{v) > (1 + 6)X{n - 1)} n {T, < t}) < e 



.1x2 
3 



<5^A(n-l) 



Now ([IT]) follows. ■ 

We shall need to upper bound the number of saturated links around any 
given node, as in the following lemma. 

Lemma 3.3 Let n and D he positive integers, and let k > 4:p£){dX){n — 1). 
Then for each t >0 

Pr({Sf (at v)>k}A{T< t}) < 2exp ( (18) 



and 



Pr({5f (ma v)>k}A{T< t}) < 2 exp ( - A ) . (19) 



Observe that if 5 > and D = D[n) —>■ 00 then for n sufficiently large we 
may for example take k as 6n. 

Proof. We use the coupling of the three processes (Xf), (Xf) and (Yt) de- 
scribed in the proof of Lemma [3.21 Let f G ^ be a node. Note that, for each 
u V, the load Yt{vu) of link vu is a Poisson random variable with mean 

A n - 2 ^ / \-, < dX. 

^ ' {n-2Y 

We write Sf{a,t v) to denote the set of links vw for calls at v that have load at 
least D at time t in the stationary superprocess (Yj); and we write Sf^ai v) = 
\S^{a.t v)\. Also, for w G W, 5° (via w) denotes the set of links uw for calls 
at some node u, and routed via w, that have load at least D at time t in 
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(Yt); and ^P(via w) = |>S°(via w)\. Then E[^f(at f)] < (n - l)pi5(rfA) and 
E[^P(via w)] < (n - l)j9D(rfA) for all times t > 0. 

For a given v E V, the loads Ff(fu) of links vu for u ^ v are determined 
by a set of {n — l)(n — 2)^^ independent Poisson random variables each with 
mean \/{n — 2^ (corresponding to — 1 choices of the other end node w and 
{n — 2Y choices of d routes for a call with end nodes v and w), and so there 
is strong concentration of measure. Note that the median m{v) of >S'P(at v) 
is at most 2{n — l)pD{dX). We can use Talagrand's inequality Lemma \2.1\ 
with c = d and r = D. This gives, for all t > 0, 

Pr(S?(at .) > ,n(v) + .) < 2exp (- ,,.^(4) + .) )- 

Now take z > 2(n — l)p£,{dX) > m{v), so that 

Pr(5f (at v) > 2z) < 2exp (-^). 

Similarly, given w eV , the loads Yt{uw) of links uw for u ^ w are determined 
by a set of ("2 ^) [(n — 2)*^ — (n — S)'^] independent random variables Po{\/ {n — 
2Y) (corresponding to calls for all possible pairs of distinct nodes v,u E 
V \ {w} choosing a route via node w). Applying Talagrand's inequality with 
c = 2 and r = D, we have, for t > and z > 2{n — l)pij{dX), 

Pr(Sf (via w) > 2z) < 2exp (-^). 

On the event {T < t} 

Xt{vu) < Xt{vu) < Yt{vu) 
for each link vu, and we deduce that inequalities f|T8|) and (|T9l) hold. ■ 

The next lemma uses the last one to show that for any node v E V, for all 
sufficiently large times t, Xt{v) is unlikely to deviate much below X{n — 1). 

Lemma 3.4 Let < 6 < 1 be a constant. Let the capacity D = D{n) he 
such that po{dX) < 5/32 and D = o{n). Let Bt be the event that Xt{v) < 
(1 — 5) A (n — 1) for some vertex v. Then there exists a constant rj = ri{6) > 
such that, for each positive integer n we have 

Vr{Bt, n {T < ti}) < 26-""/^ (20) 

for all times ti > and ^2 > ^1 + ln(4/5). 
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Proof. Observe that the left hand side of fl20p is non- decreasing in ti, and 
so we need only consider the case when t2 = ti + ln(4/(5). For each ti > 0, 
let t2 = tl + ln(4/5); and each v E V let Atj^{v) be the event that 

5° (at v) <in- 2)5/4 for all t e [h, ts] 

and let At^ = n„Ati(t'). By (fT8|) in Lemma [3731 we have 

Pr({5f (at v)>{n- 2)6/8} n {T < tj) < 2exp (-^^^) • (21) 

We now apply Lemma [2751 part (a), with a = h = {n — 2)5/8, a = 6/dX, and 
as the event that T < ti. Thus 

Pr({T < n AA^) < (i^^^^ + l)(/i(a) + 6), 

a 

where h{a) is the right hand side of ([2T1) . and 

^ = Pr(Po(A(n - l)cT > 6 + 1) < e-'5("-i)/3<i^ 
by inequality ([2]). Now summing over all w G we obtain 

ijri-2)6\ _s(n-i)/3d 



Pr({T < tijnAJ < n (rfln(4/5)A/5 + 1) (^2exp (^-i^^^ ) + 

Thus there exists a constant t] = r]{S) such that for all n G N and all ti > 
we have 

Pr({T<ti}n A7) < 2e-^"/^. 

Now, by f[TTl) in Lemma [3.11 on the event that >S'j^(at u) < {n — 2)5/4 
for each vertex u, the probability, conditional on (pt- and on the end points 
of the call, that a new call arriving at time t would not be accepted is at 
most 5/2, and so the rate at which calls for a given node v are accepted is 
at least {n — 1)A(1 — 5/2). We may now apply Lemma 12.61 in the special 
case with = 1, fi{x) as the number of calls in progress with one end v, 
and Ai = (n - 1)A(1 - 5/2). Also 1 - e'^^^"*!) = 1 - e-^'^WS) = i _ 
Hence, on Aj^, Xt^iv) stochastically dominates a Poisson random variable 
Po{{n — 1)A(1 — 35/4)). Then, using the concentration inequality ([3]), for 
all f , 

Pr(Ai n {Xt,{v) < (1 - 5)\{n - 1)}) < e-^'^("-i)/i28_ 
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Combining the above estimates we have that Pr({T < ti} fl {Xt2{v) < 
(1 — 6)\{n — 1)}) is at most 

Pr ({T < t^} n + Pr(Ai H {Xt,{v) < (1 - 6)X{n - 1)}) < 26"''"/^ 

for all n with a suitable new value of 77: thus 

Pr({r < ti} n {XM < (1 - 5)A(n - 1)}) < 26-""/^, (22) 

and the lemma follows easily. I 

To end this section we put together two of the results above to show that 
we are unlikely to observe large deviations of Xt[v) from X[n — 1) for any 
node V even during very long time intervals. 

Lemma 3.5 Given < 5 < 1, there exists a constant (3 = l3{6) > such that 
the following holds. Let the capacity D = D{n) he such that pn^dX) < 5/32, 
and say D = 0{n^). Let 6 > and let to = ln(E ||Xo||i +1) + {9 + 1) Inn. 
Let Ct denote the event that \Xt{v) — A(n — 1)| > 5\{n — 1) for some vertex v. 
Then for each positive integer n and each time ti > to + ln(4/5) 

Pr(Ct holds for some te[tiM + e''"/^]) 

Proof. Observe that Ct = AtUBt, where At and Bt are defined in Lemmas l3.2l 
and 13.41 respectively. Hence by these lemmas there exists a constant 7 > 
such that for any time t > to + ln(4/5) we have 

Pr(Ctn{T < to}) < 2e~^"/^. 

Let (3 = 7/3. We may now apply case (a) of Lemma [2.5[ with r = e^"" and 
a = 6/4, and b = 6\{n — l)/2. We use inequality ([7]) with a = (1 — 5)A(?2— 1) 
for deviations below the mean, and inequality ([H]) with a = (1 — 5/2)\{n — 1) 
for deviations above the mean. Thus for all positive integers n and all times 
ti > to + ln(4/5) we have 

Pr {{Ct for some t G [ti, ti + e^"/^]} n {T < to}) < 26"^"/^. 

We may now use f|T5l) and f|T6l) to complete the proof. ■ 
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4 Proof of Theorem 



1.1 



Recall that T denotes the departure time of the last one among those calls 
that were present in the system at time 0. Let ri = ri(n) be any time such 
that 

ri >ln(E||Xo||i +l) + (ir + 3)lnn. (23) 

4.1 Theorem HJJ: upper bound 

Since D ~ alnra/ In Inn, we have p£){d\) = n~°'~^°^^\ Then 

Pr(T > ri) < E||Xo||i e'^' < n-^-\ 

Let Na be the number of calls that arrive in the interval [ri, ri + n^]. Thus 
A^^ ~ Po{\(2)n^ ■ Let Np be the number of calls that fail in the interval 
[ti, Ti + n^]. We must show that ENp = o(l). 

Suppose first that K < d — 2 and a > {K + 2)/d. We are going to 
upper bound At in order to use ffTTj) . and for that we argue as in the proof 
of Lemma [3 -Si By Talagrand's inequality (Lemma 12. ip . for each v & V, 

Pr {{dt{v) > 4(n - l)pD{d\) + In^ n} n {T < t}) 
< Pr (Sf (at v)>4:{n- I)pd(c^A) + In^n) 
= exp{-n{\n^n/D)) = exp{-n{\n^n)). 

For t' >t > 0, let be the event that < A{n - l)pr){dX) + In^ n for all 
s G [t,t'). By the above inequality and Lemma [231 

Pr(A;;;;;;^n {T < t}) = exp{-n{\n^n)), 

and it follows using (fT5l) and (fT6l) that Pr{AT.^ ,r-i+nK) = o{n~^~^). 
By Lemma WA\ inequality (fTTj) . on A^ j/ fl {Tj = t'} 

P ^7 f -1 u w / 8(n-l)p^(rfA) + 21n'n y _^_2 
Pr(Zj fails < I — 1 = po = c>(n ). 



Let = \2ENa\ = \2\{l)n^]. Then 

ENp < E[NaIn^>No] + NoPo + E 



Npl-j I 



< o(l) + E 



TI . T-i ^-n.-'^ 



< o(l)+AroPr(A,,,,^+„A-; 
= o(l). 
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This completes the proof of the case K < d — 2 and a > {K + 2)/d. In 
the case where K > d — 2 and a > K + 3 — d, the proof is somewhat longer. 
Note first that a > 1. By Lemma [3. II inequality f|T2|) . for each time t > 0, on 



Pr {Zj fails n {T < t}\^t 



< 



n(n - 2¥ ^ 



dt-{vYl 



T<t- 



Recall from Section [3] that, for all v & V, on T < t, the number of full links 
ending in v is stochastically dominated by the number of links vu, u ^ v 
such that Yt{vu) > D, where Yf is a stationary copy of the superprocess. Let 
us consider time t = say, and call this quantity d{v). Therefore for each 
time t > Ti, on Tj = t we have 

Pr{Zj fails) < Pr(r > t) + Pr{Zj fails n {T < t}) 

2d+i 



< Pr{T>n) + 



n{n — 2 
2d+i 



-5^E[d,_(t;)%<,] 



n{n - 2y 



(24) 



Consider a fixed node v & V. We show next that 



l-Q + o(l) 



(25) 



E[(i(v)^] = E[rf>)](l + o(l)) = n 

Consider the superprocess at time 0. Let Ui,...,Ud be distinct nodes in 
V \ {v}. Let N{ui) be the number of live calls with one end v that have 
selected the link vui but none of the links vuj for j ^ i. Let N be the 
number of live calls that have selected at least two of the links vUi. Then the 
N{ui) are iid, each is Poisson with mean at most Xd, and N is Poisson with 
mean 0{l/n). 

Let X = d + a, and let A be the event that N < x. Note that Pr{A) = 
0{n~^)\ and 



E 



Yq{vu,)>D^A 



< E 



i=l 



N{ui)>D-x 



Pr{N{ui) >D-x) 



Now let Ofc be the number of partitions of 1, . . . , into exactly k non-empty 
blocks. In the sums below the Wj run over V \ {f }. We find 



E 



n I ^^Mv^j)>D^^ 
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d 

= E E[niyo(..,)>DlA] 

Wl,...,Wa j = l 

d k 

= E afe(n - l)fcE[JJ Iyo{vu,)>dU 

k=l i=l 
d 

< E[d{v)lA] + J2^kn''Pr{N{ui)>D-x)'' 

k=2 
d 

a+o{l)\k\ 



< E[d{v)] + 0{J2in' 



E[d{v)] + n 

^l-a+o(l)_ 



k=2 
-2(a-l)+o(l) 



Also 

E[d{v)%] < n'^Pr{A) = 0(n'^-^) = 0(n-"), 
and so holds, as desired. Thus on Tj > ti 

nd+l 

Pr(Zj fails ) < 0(n-^-^) + — -r n n^-"+°(^) 

n{n — 2Y 

= o(n-^-^) 
since a > K + 3 — d. It follows that 

E[Nf] < E[NAlN^>No]+o{Non-''-^) = o{l). 

as required. 

Note that a much easier proof works in the case d = 1. Let a > K + 2. 
On Tj > Ti, by symmetry 

Pr{Zj fails ) < 2Pr{Yo{vu) > D) + Pr(T > ri) 

for any pair of distinct nodes -u, f G V"; and so 

Pr{Z.j fails ) < + o{n-^-'^) = o{n-^-^). 

It now follows as above that E[A''j7] = o(l). 
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4.2 Theorem 11.11 : lower bound 



Suppose first that K < d — 2 and < a < {K + 2)/d; or that K > d — 2 and 
< a < 1. Let D ~ a Inn/ In Inn. Let < 5 < (K + 2) /d - a. We shall use 
Lemma 13.11 inequality f|T3l) to obtain a lower bound on the probability that 
a call Zf is lost; this will entail lower bounding the quantity dt-{vY- 

For < ti < let A\ ^^ be the event that 5*° (at v) < {n - 2)5/2 for all 
vertices v and all times t G [^1,^2]- Then for any v & V , any link vj and any 
time t G [^1,^2], on the probability that a call for a node f arriving at 

time t which selects link vj as its first choice is blocked by the 'partner' link 
uj (where u is the random other end of the call) is at most 5/2. 

Fix a node v. We apply Lemma with N = n — 1. For each load vector 
X and each node j 7^ f , we let fj{x) be the number of calls in progress on the 
link vj. Also, for each j we let \j = A(l — 5/2) and rij = D. It follows that 
on Aj^ ^2 the random variable 5'j^(at t>) stochastically dominates Sj^^f^^at v), 
where 

'S'f (at t;) = X]%/"^''=D' 

and the are independent immigration-death processes each with arrival 
rate X{l — 5/2), death rate 1, population at time ti and population limit D. 
It is well known that in equilibrium the n — 1 immigration-death processes 
y/""''' are iid random variables with a Poisson distribution Po{X{l — 5/2)) 
truncated at D, that is the Erlang distribution with parameters A and D. 
Since, by standard theory, each y}^''^ converges to equilibrium exponentially 
fast, there exists a constant c > such that, for all t > clnn and all j ^ v 
we have Pr(F/''-''^ = D) > n-"+°(i), and so E[^f (at v)] > n^-'^+^W _ 

We have assumed that < a < 1, and so E[S'/^(at v)] tends to infinity 
as n — s> 00. Let 0<e<K + 2 — da and refine the condition on 5 so that it 
now must satisfy 5 < — a — ^. Using inequality ([3]), on Aj_^ 

Pr(5f (at v) < n^-"-^) < exp(-n^-"+°(^)) 

for all t such that ti + clnn < t < ^2- For < ti < ^2 let ^t^.tj denote the 
event that dt{v) > n^'"'-^ for all f G and all t G [ti,t2]- 

Recall that ti was introduced in fl23l) . Let T2 = ti + clnn. Let / denote 
the interval [t2, T2 +n^], and let A denote the event Aj-^^T^^ni^ . By the above 
and Lemma 12.51 

Pr(An<,,^+„K)=o(n-^). 
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Also, we know by Lemmas 13.31 and 12.51 that 
Further by f|T6l) 

Pr(T > n) < E||Xo||i e-^' = o(n-^). 
It thus follows that 

Pr(A) = o{n-^). 

By Lemma [3.11 equation (1141) . for each t G /, on A fl {Tj = t} we have 

Pr{Zj fails | > = pq. 

Let F be the event that fewer than n'^ calls arriving during the interval / 
fail; and let Np be the number of calls that fail during this interval. Let 

Pr(F) < Pr(A) + Pr{Po{2No) < Nq) + Pr{Bin{No,po) < n') = o(n"^) 

and hence ENp = Q{n'^). 

Now consider the last case remaining, when K > d — 2 and 1 < a < 
K + 3 — d. As in the case < a < 1 considered above, for each t > T2 and 
V eV, E[dt{v)] > n^-°+°(^\ By Lemma O inequality ([13]), for each t > 0, 
on Tj = T2 + t 

Pr{Zj fails ) > n-\n - 2)-^^ ^ E[rf(,,+t)_(t;)^] 
> n-'-'J2^[d^r'^+^yiv)] 

V 

and hence EiVp = Q,(n^^^^'^^°'^°^^^) = ^1(71"), as required. 



5 Proof of Theorem 1.2 



Fix an integer d > 2 and a constant K > 0. Fix a constant < 5 < 1. 
Let 6' be a constant with 6 > (100 + K)/ In 2. We now define times tq, ti, T2 
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depending on n (not quite as when we introduced ti in (1251) in the last 
section): we let 

To > ln(E||Xo||i +1) + ^lnn, ti = tq + 9\nn smd T2 = n + . (26) 
For each t G [tq, T2\, let be the event 

{(1 - 5)X{n - 1) < Xs{v) < (1 + 5)A(n - 1) Vs G [tq, t], Vt;} ; 
by Lemma [3.51 'Py[A^^) = o{n~^~^). Also, let A] be the event 
{5f (via t;) < (n - 2)5/4 Vs G [tq, t], Vt;} ; 

then Pr{Al^) = o{n~^~^), by Lemmas 13.31 and Lemma [2. 5[ 

Recall that for each link vw, Xt{vw) is the load of link vw at time t, that 
is the number of channels in use at time t. For h = 0,1, . . ., let Lt{v, h) be the 
number of links vw at v with Xt{vw) > h (so, in particular, Lt{v, 0) = n — 1 
for all t). For an integer j, if Tj < t and the call Zj is for node v routed via 
node w and is still in progress at time t, then the height TCt{j, v) at v at time 
t of the call is one plus the number of calls in progress on vw at time t that 
arrived before it; and the call has height at v at time t if the conditions do 
not all hold. Thus we have 'Ht{j,v) < Xtiyw). For h = 1,2, . . ., we define 
Ht{v, h) to be the total number of calls on the links vw for w ^ V \ {v} with 
height at least h at time t. Clearly, Lt{v, h) < Htiy, h) for each node v and 
each positive integer h. 

Let c = max{ci,C2}, where ci and C2 are constants respectively defined 
in Sections l5.ll and [5^ below. 

5.1 Upper bound 

Let the constant ci = ci{\,d,K) be as in (l29l) below, and let D = D{n) > 
^s^s^ + Ci. We shall show that aas no calls arriving during the interval [ti, T2] 
of length fail. 

Given a positive integer Hq, for h = Hq, hi + 1, . . . let Bt{h,a) denote 
the event that Ht{v,h) < a for each v. Also, given numbers an for h = 
ho, ho + 1, . . . and times th for h = ho, ho + 1, ■ ■ ■ satisfying To <th < ti, let 

B{ho) = {Lt{v,h) < 2ah, Vt G [th,,T2], ^v}. 
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and for h = ho + 1, ho + 2, . . . let 

Bih) = {Htiv, h) < 2ah \lt e [U, T2], Vt;}. 

The idea of the proof is to choose a sequence of about In In ra/ In d numbers 
ah decreasing quickly from a constant multiple of n to zero, and an increasing 
sequence of times th for h = 0,1,2, .. . satisfying tq < < ti for all h. Then 
the aim is to show that B{ho) holds aas, and if B{h) holds aas then so does 
B{h + 1); and to deduce that B{h) holds aas for some h < D with ah = 0. 
Thus aas no link is ever saturated during [th, T2], and so no call can fail during 
that interval. 

Let = [max{8A, 768A^}] . We choose a decreasing sequence of numbers 
a/i > as follows. First, let 

. ,n — 1 n — 1 , 

ahn = mmj , k 

° ^ 8 768A ^ 

Note that aho > A(n — l)//io- Hence, on A^^, for each t G [To,r2], since 
Xt{v) < 2\{n - 1) we have Lt{v, ho) < 2\{n - l)/ho < 2aho; and so C 
B{ho). Next let 

ah f 8ah-i^ 



6A , (27) 

n — I \n — 1 J 

for h = ho + l,ho + 2, . . ., until ah < 1^{K + 2) Inn. [We shall see shortly 
that there is such an h.] When this first occurs, we let h* = h*{n) be 
the current value of h and increase ah* to 14:{K + 2) Inn. Finally we set 
ah*+i = 2K + 5. Observe that on B{h* + 1), for each t G \th*+i, T2] we have 
maXt; Xt{v) < h* + 2K + 10. Note that the recurrence fl27|) can be rewritten 
as 

= 48A-«^_i, (28) 
where oih = 8ah/ {n — 1). It follows that for hQ + l<h<h* — l 

ah = {48XY+''+-+''-"'-'af;'' < (48A ■ 

since aho — 1- -^^^ now, since 48 A ■ aho ^ ^, ioi ho < h < h* — 1 we have 

8ah ^h-ho_. 



, dh < (0.5)^^ 
n — 1 

and so h*{n) = In Inn/ Inrf + 0(1). We now set 



Ci = snp{h*{k)+AK +11-^^}, (29) 
k Ind 
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so that 

^ In In n , * , s 

D > — — + ci > h*(n) +4K + 11. 
met 

Now define an increasing sequence th of times as follows. Let 'jh = 
48("log2 (2a/j/a/i+i)] for h = ho,...,h* — 2, let jh*-! = 481og2n, and let 
7/i* = {K + 3)lian. Note that •jh*-! > 48[log2 {2ah*-i/ah*)] for n sufficiently 
large. Let t/^p = tq, and let th = th-i + ih-i ioi h = + 1, + 2, . . . , h* + 1. 
Thus th*+i = To + Yl!h=ho '^h- shall show that with high probability 
B{h* + 1) holds and so throughout the interval [t/i*+i,r2] there are no full 
links. 

Note that 

h*-2 h*-2 

th*^i-ro = ^ 7h = 48 ^ |"log2 (2aft/a/,+i)] 

h=ho h=ho 

h*-2 

< 9Q{h* - ho) + 48 ^ (log2 at - logs an+i) 

h=ho 

< 96h* + 48 log2 < 49 log2 n 
for n sufficiently large. It follows that 

th*+i - T-Q = th*-i - To + 7h._i + 7h* < 6'lnn 

for n large enough, as 6* > (100 + K) / In 2. 

Recall that C B{hQ), as we noted earlier. We shall show that 
Pr{B{h) n B{h - 1)) is small for each h = ho + 1, . . . , h* + 1, which will 
yield that Pr{B{h* + 1) is close to 1. Hence, as we discussed earlier, aas 
throughout [th*+i,T2] there are no full links. Since < ti, this shows 

that 

Pr[3j : {Tj G [n, ts]} H {Z^- is blocked}] < Pr(E(/i* + 1)) = o(n-'=). (30) 

This yields the desired upper bound of Theorem 11.21 

The main step is to prove that Pr{B{h) fl B{h — 1)) is small for each h. 
With this aim in mind, we first show that if B{h — 1) holds then aas for each 
V there exists a time th{v) G (th-i,th] such that Ht^[y){v,h) < a^. We then 
show that aas Ht{v, h) < 2ah for all t G {th{v), T2] and all v & V. For each 
node V & V and for each integer h = ho + 1, . . . , h* + 1, let 

C{v,h) = {3th{v) G {th-i,th] : Ht^^^){v,h) < ah}. 
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Let also C{h) = HyC^v, h), so that 

Ci}i) = {3w: Ht{w,h)>ah Vt G t/,]} 

is the event that there is a node u such that the number of calls with height 
at least h at u is greater than ah throughout 

Lemma 5.1 Uniformly over all h = ho + 1, . . . , h* + 1 , 

Pr(C(h)nB{h - 1)) = o{n-^-^). 

Proof. Fix a node v and a height h with Hq + I < h < h*. Let Jo{v) = th-i, 
and enumerate the jump times of the process of arrivals (possibly failing) 
and terminations of calls with one end v after time Jo('^) as Ji(f), J2{v), . . .. 
For A; = 0, 1, . . . let = Hj^^^^ (f , h) and for = 1, 2, . . . let = R^ — Rk-i, 
so that 

k 

Rk = Rq + Yj- 

i=i 

Note that each G { — 1,0,1} and is -measurable, and that the sum 
'^k:th-i<Jkiv)<th is the net change in Ht{v, h) during the interval {th-i, th]. 
Foi h = ho, . . . , h* - 1, let ruh = \12\n] \log2{2ah/ah+i)] < l^hH^ - 1) for 
n > 2. Note that for each h = Hq + 1, . . . , h* + 1 we have Jm,i_i {v) < th aas, 
since by inequality ([3]) 

Pi{Jm,_M) > th) < Pr{Po{X{n-l)jh-i) < mh-i) < e-^'^-"("-i)/8, 

which is 6"*^*^"^ 

For = 0, 1, . . . let 

Ek = n - 1) = v4°,(,) n - 1); 

and let E = n'^JlQ^'^Ek. We saw earlier that Pr(Afj = o{n-^-^). Thus 

Pr(EnB{h - 1)) < Pr{Jm,_,{v) > h) + Pr(Al) = o(n-^-i). 

Now we obtain an upper bound for the probability of positive steps and 
a lower bound for the probability of negative steps. On Ek-i fl {Rk-i > 
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ah), upper bounding Pr(Jfc(t>) is an arrival time) by 1, we obtain 



< 



< 



2 max^ Hj^(^y)_{w, h - 1) 



8ah-i 
n — 1 



n-2 

d 



d 



(for n > 3). On A^j^f^^-^_, the probability that Jk{v) is a departure time of a 
call with one end v is at least , }, — > .„/ for n > 3. It follows that, 

A(2+fl)(ra— 1) — 3A(n— 1) — ' 

for n > 3, on A\^^^_ 



3X{n-l) -3A(n-l)' 
and so, for each y > ah, on Ek_i fl {-Rfc-i = I/}, 

We note that for h < h* 

+ f 8ah^i\'^ ah _ ,^ 

Let a = and let b = ^.x^n-i) ■ Lemma [531 with r = ah and any 

value of + 1 < f < 2a/i_i, 

Pr(E n {Hj^(,){v, h) >ahWk< mh-i}\ Hj,^,){v, h) = Tq) < 26-""/^^ 

= o(n-^-i). 

It follows that, uniformly over + 1 < h < h* , 

Pr(C(/i) n B{h - 1)) = o{n-^-^). 

Now we consider /i = h* + 1. Let Jo('^) ~ and enumerate arrival 
times of calls with one end v after time th* as J2(f ), . . .. Define m^. = 

2\{n — l)7h*, which is O(nlnn), recalling that 7/^. = {K + 3) Inn. 
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Let Nt{v) denote the number of new calls for v arriving during 
with height at v at least h* + 1 on arrival. For A; = 0, 1, ... let 



Ei 



Further let E' = n'^^^'^E',^. 

Consider the call that arrives at time J'}^{v). On El^_^ it has probability 

at most pi = (^56{K + 2)n~^\nn^ of choosing a link vw with at least h* 
calls for n large enough. Further we note that, for each positive integer r, 

Fr{Bin{mh*,Pi) > r) < (m^.pi)'^ = 0((n~"'+^(lnn)°'+^)^'). 

Then, for each integer r > K + 2, 

Prm,._^^iv)>r}nBih*)) ^ 
< Pr{Bin{mh*,pi) > r) + Pr{Po{X{n - l)-fh*) > mt*) + Pr(E' n B{h*)) 

Also, on B{h*) the probability that some call with height at least h* + 1 
present at time th* survives to time th*+i is at most 28{K + 2) Inn e~^'^* = 
0(n-^~'^). It follows that 



Pr{C{h* + 1) n B{h*)) = o(n-^-') 

as required. ■ 
We now show that for each h = ho + 1, . . . , h* + 1, aas there will be no 'ex- 
cursions' that cross upwards from ah to at least 2ah, that is Ht{v, h) cannot 
exceed 2ah during (t^{h), T2] for any v E V and any h = Hq + 1, . . . , h* + 1. 

Lemma 5.2 For ho + 1 < h < h* + 1, let B\h) = B{h - 1) n C{h). Then 
uniformly over h = ho + I, ■ . ■ , h* + 1 



Pr{B{h)nB'{h)) = o{n 



-K-l\ 



Proof. The only possible start times for a crossing are arrival times dur- 
ing [th-i,T2]. Let A^o = 2X{n — 1)t2. The probability that more than Nq 
arrivals for any given node v occur during (t/i(f),T2] is o{n~^~^). We apply 
Lemma [23] with p = g^, q = g^, a = [ah\ — 1, and 
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By Lemma [2.41 the probability that any given excursion leads to a 'crossing' 
is at most {ql / q^)^""^^'^ ■, and so for /i = /iq + 1, • • • , /i* + 1 

Pr(5(/i) n < A^n(g+/g^ + Pr(E) + oin'^-^). 

For h = ho + l,...,h* 

for /i = /i* + 1 

/C)"'"' < (3A/2ir + 5)'^+=^(112(i^ + 2) lnn)'^(2i^+3)^^ _ ^)-{d-i)(2A'+3) 
= 0(n-2^-='(lnn)^^+'^), 

and the lemma follows. ■ 

We may now complete the proof. As C Bh^ and Pr{A'^^) = o{n~^~^), 
we have 

Pr{B{h* + 1)) 

< Pr(I^) + Pr(fi()^nA°J+ ^ Pr(5()^n5(/i-l)) 

h=ho+l 

h*+l h*+l 

= Pr(A^)+ J2 Pr(C(h) n B{h - 1)) + J2 Pr(B(l^nC{h)nB{h- 

h=ho+l h=ho+l 

= o(n-^). 

This completes the proof of (1501) and thus of the upper bound of Theorem 11.21 
5.2 Lower bound 

Let < e < min{l, (i^ + l)/d}. Recall that 6' is a constant satisfying 
9 > (100 + _ft')/ln2. Let the constant C2 = C2{\,d,K) be as defined below, 
and let D = D{n) < — C2. Recall that tq, ti and T2 are defined in ([2HD 
at the start of this section. We consider the interval [Ti,r2] of length . 
We shall show that aas for each v at least (n — l)^"*^ links vw incident on v 
are saturated (and so unavailable) throughout the interval, and hence aas at 
least n^+'^-<^-°W calls arriving during the interval fail. 

Given a non-negative integer h and a real value a > 0, let Bt{h, a) denote 
the event that Lt{v,h) > a for each v. Given also a sequence of numbers 
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ah, h = 0, 1, 2, . . . and a sequence of times th, h = 0, 1, . . ., let B{h) = 
^t£[th,T2]Bt{h, ah)', thus B{h) is the event that for every v E V, throughout 
the interval [th, T2], the number Lt{v, h) of links with one end v that carry at 
least h calls remains at least ah- We shall choose positive numbers ao,ai, . . ., 
starting with ao = n — 1 and decreasing rapidly. These will satisfy 2ah+i < 
a;/i(l — e^-*^), so that {ah — 2ah+iY — {(^h/^Y- We shall further choose an 
increasing sequence of times th, h = . . ., such that tq < th < ti for each 
h. We want to show that B{D) occurs aas, with a value ao = fl{n), so 
that there are always many saturated links. Analogously to the upper bound 
proof in Section 15. the main task is to show that Pr{B{h) fl B{h — 1)) is 
small for each h < D. 

Again similarly to the upper bound proof, we first show that on Btf^^^^t^ 
1, ah-i) aas there exists a time th{v) G {th-i, th\ such that Lf^(t,)(f , h) > 2ah- 
We then show that on B{h — 1) aas Lt{v,h) never falls below ah dur- 
ing {th{v),TQ\. The numbers ah are given as follows. We let ao = ^ ~ 1) 
and for /i = 1, 2, . . . 



ah min{l, A} / ah 



n-l 2Ae'^h \n-l 



(31) 



Thus 2ah < a/i-i(l — as required, since < e — 1. We need to choose 
the constant C2 in the upper bound on D{n) above such that for n sufficiently 
large 

It is easy to check that such a choice is possible. To see this, let u = ""^^^^-^ 
and let = Then < i/ < 1, /3o = 1 and for /i = 1, 2, . . . 

f^h = (32) 
It follows that for each positive integer h 



To upper bound the denominator, note that 

h 



ln{h{h - lY{h - 2Y' . . . 2^^'^') =d''J2 In i < c^d'' 



i=2 
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for some constant C3 > 0, and so the denominator is at most e"^^ . It follows 
that for each h & N 

Let C4 be such that ci~'^*(ln(-) + C3) < e. Hence ii h < lnln(n — l)/lnci — C4 
then 

(5n > exp(-(ln(n-l))ci-^^(ln(l/z/)+C3)) 
> exp(— eln(n — 1)) = (n — 

For h — 0, 1, . . . we let 7/, = n-TTTrT-rr- Now define an increasing 

' ' max{l,A}(/i+l) ° 

sequence of times th as follows. Let to — Tq, and for /i = 1, . . . , let = 
th-i + Ih-i- Then 

D-l D ^ 

to - To = <^'^-f^ < 41n(L' + 1) = O(lnlnlnn). 

h=0 h=l 

It follows that to < Ti for n sufficiently large. 
We now recall that 

B{h)^{Lt{v,h)>ah yte[th,T2],yv}, h=0,l,.... 

Thus Pr(£)(0)) = 1; and we prove by induction that Pr{B{h)) — o{n~^''^) 
for h = 1, . . . , D, so that aas throughout [ti, for each v there are at least 
(n — 1)^^^ saturated links vw incident on v. 

Fix a node v and an integer h > 1. Let Jo{v) = th-i and enumer- 
ate the jump times of the process after time Jo{v) that concern node v as 
Ji{v), J2(v)- ■ ■ •• For k — 0,1, ... let Rk — Lj^(^^^{v, h) and for A; = 1, 2, . . . let 
^k — Rk~ Rk-i-i SO that 

k 

Rk — Ro + ■ 

Then each Yk e {-1, 0, 1}, is 0jj^(^)-measurable, and ^k:tn-i<Jkiv)<tn is the 
net change in Lt{v, h) during {th_i,th\. Let ruh = 2min{l, A}(n — = 
\\{n— 1)7/1. Note that Jmh-A'^) ^ ^/i with high probability for h < D, since 

PriJrn,_Av) > h) < Pr(Po(A(n-l)7/.-i < m,_i) < e-^''-^("-^)/« = o{n~''-') 
For A; = 0, 1, . . . let 
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Let E = n'^^^'~'^Ek. Recall that Pr{A°^ U AjJ = o(n-^-i). Thus 



Pr{E n - 1)) < Pr{J^,_,{v) > t,) + Pr{Al U A^, 



o(n 



Now we want to give a lower bound on the probability that Yk = 1. First 
note that on A^j^f^,^^_ the probability that Jk{v) is an arrival time (for a call 
for v) is at least 1/(2 + 5) > |. Now consider picking the random other 
end u of the call and the random intermediate nodes Wi, . . . , in the order 
wi, u, W2, ■ ■ ■ , Wd- On we have S'j^(^)_(via w) < (n — 2)/2 for all nodes 

w; and so, whatever wi is picked, the probability that uwi is saturated is at 
most |. Hence on fl A^j^^^-^_ we have 

Pr(n = l|0j,(,)_) 

> 1 ^JfcM-(^>^-l)-^Jfc(f)-(^^>^) 1 /^^JfcW-(^>^-l)-l-^JfcW-(^^,^) 



> 



3 n-1 2 \ n-2 

1 / Lj^^^)_{v,h- 1) - 1 - Lj^^^y{v,h)Y 



d-l 



6 V n-1 



Similarly on A^j^^^-^_ the probability that Jk{v) is a departure time of a given 
call with one end v is at most w„ J:, — tt < tt-^— rr, and so 

A(2— fl)(n— 1) — A(n— 1) ' 

pr(n = < "^^""-'"'"^(^^y-'"-'^^" . 

It follows that on Ek-i fl {Rk~i < 2ah) 

Also for each y < 2ah, on E^-i fl (i?A;_i = |/) 

T> /V iiA \ ^ hy 2hah 



A(n- 1) - A(n- 1) 
^ 2/iafe _ _ 

- min{l, A}(n - 1) ~ ' 

We note that qj^ = |g,| for each positive integer h. 
For each node G V and each positive integer h let 

C{v,h) = {3th{v) e (t/i-i,i/i] : Lt^^v){v,h) > 2ah}, 
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and let C{h) = r\yC{v, h). We now show that, uniformly over h = 1, . . . , D, 

Pr(C(h)nB{h - 1)) = o{n-^-^). 

Let p = ql, let ri = 2ah, and let tq be any positive integer less than 2ah- 
Note that q^nih-i > 4q;/i > 2(ri — ro). By a natural 'reversed' version of 
Lemma 12.21 

Pr{E n (Lj^(^) (t;, /i) < 2a/, VA; G {1, , . . . , m/,_i) | Ljo(^) (t;, h - 1) = ro) 
It follows that 

Pr(C{fi}nB{h - 1)) = o{n-^-^). 

We now need to prove that for each h = 1,2, . . . , D, aas there will be no 
excursions that cross downwards from 2ah to less than ah, that is none of 
the numbers Lt{v, h) can drop below au during {t^{h), T2\. 

For each positive integer h let B'{h) = B{h — 1) fl C{h). We shall show 
that uniformly over 1 < h < D 

Pr(B(h) n B'{h)) = o{n-^-^). 

The only possible start times for a crossing are departure times affecting 
links at v during [th-i,T2]. Let Nq = 4Ar2(n — 1). Let F denote the event 
that there are more than A^^o such departures in [t/i_i,T2]. Then 

Pr(FnA°J = o(n-^-i). 

We apply a reversed version of Lemma [XH with p = qj^ , q = , = [cih\ — 1, 
and 

E, = n n - !)• 

The probability that any given excursion leads to a 'crossing' is at most 
iQh/Qh)^'"'^'^ < (l/2)"'>-2. It follows that 

Pr{B{h)nB'{h)) < nA/'o(0.5)"'^~2 + nPr(Fn + Pr(E) 
< AXT2{n - l)n(0.5)"'^~2 + o{n-^-^) 
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The proof may now be completed in the same way as the proof of the upper 
bound. We have 



D 



Pr{B{D)) < Pr(fi(0)) + ^ Pr{B{h) n B{h - 1)) 

h=l 

D D 

= P^cJh) n B{h + P^Wh) n c{h) n B{h - 1)) 



h=l h=l 

= o(n-^). 

For t G [ti, let Ft be the event that for each v at least [n — links vw 
with one end v stay saturated throughout the interval [Ti,t]. By the above, 
provided that D < ^^j^ — C2 the event Fj.^ holds aas. Now let s e [Ti,r2] 
and consider a new call arriving at time s. On Fg_ the probability that this 
call is blocked is at least 

Pi = > -ra 

for n sufficiently large. Let Ni = l^X^^^^n^]. Let B* be the event that fewer 
than b* = ^An^+^~'^^ calls are lost. Then 

Pr{B*) < Pr(F^^)+Pr{Po{x(^n^) < Ni) +Pr{Bin{Ni,pi) < b*) = o(l). 
This completes the proof of the lower bound of Theorem 11.21 



6 Concluding remarks 

We have considered the performance of two algorithms for a continuous-time 
network routing problem, strengthening and extending the earlier results 
in [9] and [1]. 

The analysis in [9] (see also [Ij) suggests that the performance of the 
model can be upper and lower bounded by differential equations. While this 
analysis is non-rigorous, it is hoped that a suitable differential equation ap- 
proximation, and concentration of measure bounds, can indeed be obtained. 
The main challenge is to disentangle the complex dependencies within sub- 
sets of links to obtain a tractable asymptotic approximation for the generator 
of the underlying Markov process. The details will appear in [5]. 
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For simplicity we have assumed throughout that the underlying network 
is a complete graph, but our results will carry over in a straightforward way 
to a suitably 'dense' subnetwork. Consider for example the upper bound 
in Theorem 11.21 part (a). Let 6 > 0, and suppose that, in the network 
with n vertices, for each pair of distinct vertices u and v the number of 
possible intermediate nodes is at least 6n. Then minor alterations to the 
proof of Theorem 11.21 part (a) show that we obtain the same conclusion: if 
D{n) > Inlnn/lnci + c and we use the BDAR algorithm, then the expected 
number of failing calls during the interval of length is o(l). The only 
difference is that now the constant c depends also on 6. Note that the leading 
term Inlnn/lnc? depends only on the problem size n and the number d of 
choices, and not on 6 (or on A or K). 
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